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We review current best estimates of the strength and detectability of the gravi- 
tational waves from a variety of sources, for both ground-based and space-based 
detectors, and we describe the information carried by the waves. 

1 Preamble 

Major relativity conferences such as GR16 traditionally include an overview 
talk on gravitationaLwave (GLW) sources. Some excellent recent ones include 
those by Flanagan El, Finn □, and Bender et al. Q. Such talks are always 
by theorists, and can be described as basically informed speculation, since 
until now detectors sufficiently sensitive to detect GW's have not existed. 
But km-scale laser interferometers (IFO's) are now coming on-line, and it 
seems very likely these will detect mergers of compact binaries within the 
next seven years, and possibly much sooner. There are several other classes of 
sources that the ground-based IFO's might well detect: massive star collapse 
(supernovae and hypernovae), rapidly rotating neutron stars, and possibly 
a GW stochastic background created in the early universe. A space-based 
interferometer, LISA, is also planned (though not yet fully funded), and could 
fly in ~ 2011. There is one type of source — short-period galactic binaries 
— that LISA is guaranteed to observe (at its planned sensitivity), plus a list 
of very promising candidates: the inspiral and merger of supermassive black 
holes (SMBH's), the inspiral and capture of compact objects by SMBH's, and 
sources in the very early universe. 

In this article we review the various GW sources that have been studied, 
for both ground and space-based detectors, summarizing the best available 
estimates for rates and source strengths and describing some of the informa- 
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tion that may be extracted from the waves. We wiU try to emphasize what is 
new, but much here is not, and indeed many sections of this papeii-have been 
adapted, with minor updating, from an earher review by Thorne □ or from a 
pubHcly available LIGO-Project document by Thorne □ (based on extensive 
input from the colleagues listed as consultants at the beginning of this article). 

Gravitational-wave detection efforts focus on four frequency bands: (i) 
The extremely low frequency (ELF) band, 10^^^ to 10"^* Hz, in which waves 
will be sought via their impriril on the polarization of the cosmic microwave 
background (CMB) radiation and the only expected source is primordial 
gravitational fluctuations amplified by inflation of the universe, (ii) The very 
low frequency (VLF) band, 10~^ to 10~^ Hz, in which the waves are sought 
using high-stability pulsar timing 13, and the expected sources are processes 
in the very early universe and extremely massive binary black hole systems, 
(iii) The low frequency (LF) band, 10""* to 1 Hz, in which LISA will operate, 
and (iv) the high frequency (HF) band, 1 to 10^ Hz, in which earth-based 
detectors (interferometers and resonant-mass detectors) operate. 

In this article we shall confine attention to the LF and HF bands, since 
that is where the richest variety of sources reside, that is where the vast 
majority of the current experimental activity concentrates, and that is where 
the first firm detections are likely to be made. We shall not say much about 
details or current status or plans for gravitational wave detection, as there 
are matty other reviess/ articles devoted to this, e.g. those by Kawamura 13, 
Schutz El, Robertson and Barish and Weiss ll^. 

In Sec. H we shall discuss high-frequency sources being sought by earth- 
based detectors, and in Sec. ^, low-frequency sources that will be sought by 
LISA (ancLspme of which are already being sought by Doppler tracking of 
spacecraft t3) . 

2 High-Frequency Sources Sought by Earth-Based Detectors 

2.1 Fiducial Noise Curves 

In discussing gravitational-wave sources, it is important to put them in the 
context of detector sensitivities, both near-term and planned. As is discussed 
in Kawamura's chapter in this volume El, the first-generation kilometer-scale 
interferometric detectors (LIGO, VIRGO, GEO600, TAMA300) were being 
assembled at the time of GR16 and will begin their initial GW searches about 
the time that our review is published. In about 2007 these first-generation 
interferometers (IFOs) will be upgraded to produce second-generation IFOs, 
with noise curves that can be adjusted by adjusting the position and refiec- 
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tivity of one mirror. In parallel with this IFO experimental effort, a netwpijk 
of narrow-band resonant-mass detectors is operating and being improved 113. 

In this review we shall find it convenient to evaluate prospects for detec- 
tion with respect to three representative, fiducial noise curves: (i) that for 
the first-generation LIGO IFOs (the LIGO-I noise curve) ; (ii) that for the 
second-generation LIGO IFOs (sometimes called advanced LIGO IFOs and 
sometimes called LIGO-II IFOs), when they are operated in a "wide-band" 
mode (the WB LIGO-IInoise curve); and (iii) that for second-generation IFOs 
operated in a "narrow-band" mode that is optimized for sources near 600 Hz 
frequency (the NB LIGO-IInoise curve). Though focusing on LIGO smacks of 
American egocentrism, the reason is pragmatic: the first-generation VIRGO 
noise curve would largely lie between ithe initial LIGO (LIGO-I) curve and 
the advanced LIGO (LIGO-II) curves anyway, VIRGO does not publish an 
official second-generation noise curve, and GEO600 and TAMA300 are purely 
first-generation instruments. 

Our three fiducial noise curves are shown in Fig. |^ along with the wave 
strengths for a number of HF sources to be discussed below. The horizontal 
axis in the figure is GW frequency /; the vertical axis is the square root of 
the spectral density of an IFO's arm-length difference, i.e. the IFO's "strain 
per root Hertz" h{f) = •\/»S'al/l(/)- The (dimensionless) rms noise in a 
bandwidth A/ around frequency / is given by /irms — h{f)\/Af. 

Notice that, compared to LIGO I, the WB LIGO-II IFOs (i) will lower the 
amplitude noise by a factor ~ 15 at the frequencies of best sensitivity / ~ 100- 
200 Hz, and (ii) will widen the band of high sensitivity at both low frequencies 
(pushing it down to ^ 20 Hz) and high frequencies (pushing it up to ^ 1000 
Hz) . The lowered WB noise at optimal frequencies will increase the event rate 
for distant, extragalactic sources by a factor ~ IS'^ ~ 3000. Opening up lower 
and higher frequency bands will bring us into the domains of new sources: 
colliding, massive black holes and stochastic background at low frequencies; 
low-mass X-ray binaries, fast pulsars and tidal disruption of neutron stars by 
black holes at high frequencies. Noise-curve reshaping, as illustrated by the 
NB LIGO-II IFO, will reduce the noise by a factor ^3 to 5 relative to the 
WB IFO within some chosen narrow frequency band A/// ~ 0.2 in which 
targeted periodic sources (e.g. low-mass X-ray binaries) are expected to lie. 
We shall refer to this as "narrow-band tuning" of the IFO. 

The expected LIGO II noise curves depend somewhat on whether the 
mirrors' substrate is made from sapphire (the preferred material) or from 
fused silica (the backup). We assume sapphire in this article. If problems 
with sapphire force fused silica to be used, the WB IFOs will likely be a few 
tens of percent more noisy than Fig. |l| near the noise-curve minimum, and 
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Figure 1. The noise h{f) in several planned LIGO interferometers plotted as a function 
of gravity-wave frequency /, and compared with the estimated signal strengths hs{f) from 
various sources. The signal strength hs{f) is defined in such a way that, wherever a signal 
point or curve lies above the interferometer's noise curve, the signal, coming from a ran- 
dom direction on the sky and with a random orientation, is detectable with a false alarm 
probability of less than one per cent; see the text for greater detail and discussion. 



the overall GW event rate may be worse by about a factor 2. 

Figure 1 shows, along with the noise curves, the estimated signal strengths 
hs{f) for various sources. These signal strengths are defined in such a wayllj 
that the ratio hs{f)/h{f) is equal to the ratio of signal S to noise thresh- 
old T, rms averaged over source directions and orientations, hs{f)/h{f) = 
(S^/T^)^/^, with the threshold being that at which the false alarm probability 
is one per cent when using the best currently known, practical data analysis 



cutlerThorne: submitted to World Scientific on February 5, 2008 4 



algorithm. (For broad-band sources, the algorithm is assumed to integrate 
over a bandwidth equal to frequency and to use the output from LIGO's two 
4km IFOs and one 2km IPO thereby removing all non-Gaussian noise. 
For periodic sources such as spinning neutron stars, the algorithm uses data 
from only one 4 km IFO, usually narrow banded, the noise again is assumed 
Gaussian, and the signal is integrated for 10^ sec, except in cases such as Low 
Mass X-Ray Binaries where there is little gain from integrating so long.) This 
definition of hs{f) means that, wherever a signal point or signal curve lies 
above the IFO noise curve, the signal, coming from a random direction on the 
sky and with a random orientation, is detectable with a false alarm probability 
of less than one per cent. 

2.2 Overview of High- Frequency Sources 

In this section, we shall discuss the high-frequency sources briefly, in turn, and 
then in subsequent sections we shall discuss them in greater detail, focusing on 
the likelihood of detection and the science we expect to extract from detected 
waves. 

The three arrowed, long-dashed lines in Fig. 1 represent the signal hs{f) 
from neutron-star (NS) and black-hole (BH) binaries in the last few 
minutes of their inspiral, assuming masses M — 1.4Mq for each NS and 
M — IOMq for eack BH. These sources are best searched for by the method 
of matched filters E3. Using matched filters, LIGO's initial IFOs can detect 
NS /NS inspirals (with a 1 per cent false alarm probability) out to a distance 
of 20 Mpc [top arrowed line] ; the wide-band advanced (LIGO-II) IFOs can do 
so 15 times farther, out to 300 Mpc for NS/NS and out to 650 Mpc for NS/BH 
[bottom arrowed line]. For BH/BH inspiral, the wide-band IFOs can see so 
far that cosmological effects are important. For definiteness, throughout this 
proposal we assume a Hubble expansion rate Hg = 65 km/s/Mpc, a cold- 
matter density 0.4 of that required to close the universe ^Im = 0.4, and 
a vacuum energy density (cosmological constant) 0.6 of closure J^a — 0.6. 
Then the wide-band IFOs can see {IOMq /IOMq) BH/BH inspirals out to a 
cosmological redshift z — 0.4. The binary inspiral rates at these advanced IFO 
distances are likely to be many per year; see Box 1 at the end of this section. 
The middle arrowed line is the signal from BH/BH inspiral at 400Mpc, where 
the geometric mean of the event-rate estimates for BH/BH field binaries is 
three per year (third column of table in Box 1). 

The tidal disruption of a NS by its BH companion at the endpoint 
of NS/BH inspiral should produce gravitational waves that carry detailed 
information about the NS structure and equation of state. The advanced IFOs 



cutlerThorne: submitted to World Scientific on February 5, 2008 5 



may detect these waves and extract their information; see Sec. p,.^ . This tidal 
disruption is a promising candidate for the trigger of some gamma-ray 
bursts, as is the final merger of the two NS's in a NS/NS binary. A gamma- 
burst / gravitational-wave coincidence would be of great value in revealing 
the nature of gamma-burst sources; see Sec. ^.8| 

For BH/BH binaries much heavier than ~ IOA/q/IOA/o, most of the 
gravitational signal is likely to come from the black holes' merger and 
the vibrational ringdovifn of the final black hole, rather than from the 
inspiral. Rough estimates discussed in Sec. 2.5 suggest that, if the holes 
are rapidly spinning (within a few per cent of the fastest spin allowed, i.e. 
a/m ^ 0.98), then the wide-band IFOs can see the merger waves from two 
2OM0 holes out to z = 1 and two 50Mq holes out to z = 2; see the down- 
sloping, non-arrowed, dashed lines in Fig. 1. The event rates at these distances 
could well be many per year, and the waves from such mergers will carry rich 



physical and astrophysical information; see Sec. 2.5 



The triangles, star, large dots, and up-sloaping short-dashed lines in Fig. 
1 represent signals from slightly deformed, spinning neutron stars. The 
most interesting of these is a class of objects called low-mass X-ray bina- 
ries (LMXB's). These are neutron stars that are being torqued by accretion 
from a companion, but that seem to be locked into spin periods in the range 
~ 300 - 600 revolutions per second. The most plausible explanation for this 
apparent locking is that the accretion is producing an asymmetry that radi- 
ates gravitational waves, y4f]^ torque the star's spin down at the same rate 
as accretion torques it up Il3^ta. Assuming this to be true, one can deduce an 
LMXB's wave strength hs from its measured X-ray flux and its spin frequency 
[with the frequency inferred, sometimes to within a few Hz but not better, 
from nearly coherent oscillations (NCOs) in type-I X-ray bursts, or less rjeli- 
ably from frequency splittings of quasiperiodic (QPO) X-ray oscillations t2l.] 
The spin frequency, and thence the gravitational-wave frequency, will wan- 
der somewhat due to fluctuations in accretion (which can be estimated by 
monitoring the X-ray flux) and due to poorly known orbital parameters. As 
a result, in searching for an LMXB's waves one can only perform coherent 
integrations for about 20 days; thereafter, one must stack the-signals inco- 
herently, allowing for unknown shifts of the wave frequency E30. When one 
uses this "stack-slide" method of data analysis, the resulting signal strengths 
may improve only slightly for integration times longer than 20 days E3. (The 
most sensitive LMXB searches will probably be with NB IFO's, with the op- 
timum sensitivity band re-adjusted every ~ 20 days to cover the whole range 
- 400- 1000 Hz.) 

Assuming 20 days of integration using a single 4- km IFO, the estimated 
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signal strengths and frequencies for the strongest known LMXB's are shown 
by the big dots and the star in Fig. 1. The estimated strengths assume a 
steady-state balance of accretion torque by gravitational-wave torque, and 
that the GW's are generated by density inhomogneneities (as opcosed to NS 
wobble or unstable modes), so few is twice the spin frequency liZl. For every 
LMXB in Fig. 1 except the weakest one, the estimated rotation frequency 
is based on QPO splittings rather than NCOs, which means that about half 
of the frequencies might be double these estimates: ~ 1200 Hz for density 
inhomogeneities and ~ 800 Hz for r-mode sloshing, rather than ^ 600 Hz. 
Doubling the frequency above 600 Hz reduces the emission amplitude by a 
factor 2 (at fixed X-ray flux assuming a steady-state torque balance), and 

1/2 y — 

increases the amplitude noise 5*^ in an advanced NB IFO by a factor ~ v2. 
As a result, at ^ 1200 Hz frequency, Sco X-1 would still be readily detectable 
in a NB IFO but not in a WB IFO, and the strongest of the other LMXB's 
would be marginal. By contrast, if the frequencies are ~ 600 Hz or ~ 400 
Hz and the waves are in a steady state, then Sco X-1 should be very easily 
detected by a WB IFO, and several LMXB's should be detectable by narrow 
banding (NB). 

The advanced IFOs can perform interesting searches for waves from 
known radio pulsars, such as the Crab and Vela for which the current 
upper limits (based on the pulsars' observed spindown rates) are shown as 
triangles. A Crab search, using coherent integrations based on the star's ob- 
served (slightly wandering) rotation period, could improve the limit on the 
Crab's wave amplitude by a factor 100 and constrain the star's gravitational 
ellipticity to e <J 7 x 10~^ — which is approaching the realm of physically 
plausible ellipticities, e <^ 1Q~^. (We shall discuss the physical origins of 
neutron-star ellipticities and their plausible ranges in Sec. 2.7.) 

More interesting will be searches for waves from known, fast pulsars, 
since the signal strength scales as hg cx ef'^/r (where r is distance to the 
source). The up-sloping short-dashed lines in Fig. 1 show some examples of 
signal strengths. With a narrow-band IFO tuned to the vicinity of such a fast 
pulsar, the waves would be detectable when e ^ 2x 10^^(100Hz//)^(r/10kpc), 
which is in the realm of plausible ellipticities for pulsars throughout our galaxy 
so long / exceeds 200 Hz, i.e. the spin frequency exceeds 100 Hz. 

Also of great interest will be searches for previously unknown spin- 
ning neutron stars, for which the signal strengths hg will be reduced by 
a factor of a few to 15 by the lack of prior information about the fre- 
quency and its evolution and the direction to the source (which_d£termines 
the time-evolving doppler shift produced by the earth's motion) EjEj. A tun- 
able, narrow-band IFO will be crucial to such searches. One can search more 
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deeply by a factor of several, using a narrow-banded IFO that dwells on a 
given frequency and its neighborhood for a few days or weeks and then moves 
on to another frequency, than using a broad-interferometer that collects signal 
at all frequencies simultaneously for a year. Such searches will be in the band 
of physically plausible ellipticities, for stars throughout our galaxy, if / ^ 400 
Hz (spin frequency above ~ 200 revolutions per second). 

One can search for a stochastic background of gtavitational waves by 
cross correlating the outputs of LIGO's two 4 km IFOs cB. For such a search 
the signal strengths hs{f) are shown in Fig. 1 as downward-sloping dotted 
lines, assuming a cross-correlation of 4 months of (not necessarily contiguous) 
data, and isotropic waves. The lines are labeled by the waves' energy density 
n(/) in a bandwidth equal to frequency and in units of the density to close 
the universe, fl{f) = {f cLEq^^ / df) / pdosurc- Unfortunately, the frequency of 
optimal a priori sensitivity, / ~ 70 Hz (where the noise curve is parallel 
to the dotted Q, lines), is near the center of a dead band for LIGO. This 
dead band arises from the fact that l/70Hz ~ 14 ms ia_about the round-trip 
gravity-wave travel time between the two LIGO sites Efl. The result is a net 
debilitation of the stochastic background sensitivity by a factor of a few: the 
initial IFOs can detect an isotropic background with down to ^ 10~^, while 
advanced, wide-band IFOs can reach down to 17 ~ 5 x 10~^, and a factor ^ 2 
lower than this if it is reoptimized for low frequencies. These are interesting 
sensitivities, able to test a wide range of speculations about the physics of the 
early universe and perhaps even detect waves — but not able to reach the 
level predicted by standard inflation, Q. <^ 10~^^. Such a detection could have 
profound implications for physics and cosmology; see Sec. 

Other waves that advanced IFOs in LIGO will seek and may detect are 
those from the stellar core collapse that triggers supernovae and the 
boiling of the nascent neutron star at the endpoint of that collapse (Sec. 
2.6), accretion induced collapse of white dwarfs (Sec. |2.6| ), and totally 
unknown sources. 

We now give a more detailed discussion of each of the high-frequency 
sources, focusing especially on event rate estimates and the information that 
the waves may bring. Box 1 gives a brief summary of the various sources and 
their detectabihty by LIGO I and II (WB and NB). 



2.3 Inspiraling NS and BH Binaries with Mbh WMq 

As we have discussed, wide-band advanced-LIGO (LIGO-II) IFOs can detect 
the waves from NS/NS inspirals out to SOOMpc, NS/BH out to 650 Mpc, and 
BH/BH out to z — 0.4. The event rates out to these distances can be esti- 
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Box 1 

Brief Summary of Detection Capabilities of Advanced LIGO (LIGO-II) Interferometers 

• Inspiral of NS/NS, NS/BH and BH/BH Binaries: The table below 00 shows estimated rates 
TZged in our galaxy (with masses ~ IAMq for NS and ^ lOM© for BH), the distances Vi and I?wb to 
which initial IFOs and advanced WB IFOs can detect them, and corresponding estimates of detection 



rates TZi and TZwb', Sec. 2.3 



Di 

Til, yr 
-DwB 
T^WB, yr 



-1 



-1 



NS/NS 

10-6-5 X 10-4 
20 Mpc 
3 X 10-4 - 0.3 
300 Mpc 
1 - 800 



NS/BH 

<, 10-7-10 
43 Mpc 
^ 4 X 10-4 
650 Mpc 
^ 1 - 1500 



BH/BH in field BH/BH in clusters 



-4 



0.6 



<. 10-'^ 

100 

< 4 X 10 
z = 0.4 
^ 30 - 4000 



10' 



-3 



- 0.6 



- 10-6-10-5 
100 

~ 0.04 - 0.6 
z = 0.4 

- 300 - 4000 



• Tidal disruption of NS by BH in NS/BH binaries: First crude estimates suggest WB IFOs can 
measure onset of disruption at 140Mpc well enough to deduce the NS radius to 15% accuracy (compared 
to current uncertainties of a factor ^ 2); see table above for rates; Sec. 2.4. 



• BH/BH merger and ringdown: Rough estimates suggest detectability, by WB IFOs out to the 
cosmological distances shown in Fig. 2(b); rates for BH/BH total mass ~ 2OM0 are in table above; rates 
for much larger masses are unknown; Sec. 2.5. 

• Low-Mass X-Ray Binaries: If accretion's spin-up torque on NS due is counterbalanced by 
gravitational-wave-emission torque, then WB IFOs can detect Sco X-1, and NB IFOs can detect ~ 6 
other known LMXB's; Sec. 

• Fast, Known Spinning NS's (Pulsars with pulse frequency above 100 Hz): Detectable by a 
advanced NB IFO in 3 months' integration time, if NS ellipticity is e > 2 x 10-'*(1000Hz//)2(r/10kpc), 
where / is gravity wave frequency (twice the pulsar frequency) and r is distance; actual ellipticities are 
unknown, but plausible range is e lO^^; Sec. 2.7 . 



• Fast, Unknown Spinning NS's: Unknown frequency wandering and doppler shifts degrade the 
detectable ellipticity e by a factor of a few to ^ 15, so detection with a NB IFO requires e ;> (0.6 to 
3) X 10-5(100Hz//)2(r/10kpc); Sees. ^ and [I?. 



• Centrifugally Hung-Up Proto Neutron Stars in White-Dwarf Accretion-Induced Collapse 
and in Supernovae: Dynamics of star very poorly understood; if instability deforms star into tumbling 
bar, may be detectable by WB IFOs to ~ 20 Mpc (the Virgo Cluster), and possibly farther; event rates 
uncertain but could be enough for detection; Sec. 2.6. 

• Convection of Supernova Core: May be detectable by WB IFOs, via correlations with neutrinos, 
for supernovae in our Galaxy and possibly Magellanic Clouds; Sec. |2.6| . 

• Gamma Ray Bursts: If triggered by NS/BH mergers, a few per year could be detectable by WB 
IFOs; if none are seen individually, statistical studies could nevertheless confirm gravity-wave emission 
by the gamma-burst triggers; Sec. 2.8. 

• Stochastic Background: Detectable by cross correlating Hanford and Livingston 4km detector 
outputs, if = (gravitational- wave energy in A/ / ~ 40 Hz) 5 x 10^^; there are many possible 
sources of such waves in very early universe, all very speculative; Sec. |2.9|. 
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mated from observational data in our own galaxy e3, and an extrapolation out 
through the universe based on the dpisi±}2_pf massive stars (which can be de- 
duced by several different methods) Ej£j-EZI. The resulting rates (as compiled 
and evaluated by V. Kalogera Ej) are quite uncertain but very promising; see 
Box 1. 

For NS/NS binaries, the event rate in our galaxy is constrained by the 
results of radio astronomers' searches for binary pulsars that will merge, due 
to gravitational radiation emission, in_Iess than the age of the universe, and 
by other aspects of pulsar searches c3. The resulting constraints, 10~^/yr 
<, T^ged < 5 X 10^^/yr, extrapolate to a NS/NS event rate for advanced IFOs 
between 1/yr and 800/yr. Searches in our galaxy for NS/BH binaries in 
which the NS is a pulsar have failed to find any as yet, so we must turn to 
much less reliable estimates based on "population synthesis" (simulations of 
the evolution of a population of progenitor binary systems to determine the 
number that make NS /BH binaries compact enough to merge in less than the 
universe's age). Population synthesis gives aJNS/BH event rate in our galaxy 
in the range ~ 10~''/yr ^ 7?.gai <, 10~''/yr E3, though it is possible the rate 
could be even less than this (hence the <^ 10~'^/yr in Box 1). Extrapolating 
out into the universe, Kalogera finds a NS/BH rate in WB advanced IFOs 
between ^ 1/yr and about 1500/yr. A similar analysis for BH/BH binaries 
based on population synthesis (third ix)lumn of the table in Box 1) gives a 
rate between <^ 30/yr and ~ 4000 / yr Ej. That population synthesis estimate 
ignores the likely role of globular clusters aad other types of dense star clusters 
as "machines" for making BH/BH binaries single black holes, being heavier 
than most stars in a globular, sink to the center via tidal friction, find each 
other, and make binaries; then the BH/BH binaries get "hardened" (made 
more compact) by interaction with other black holes, reaching sizes where 
gravitational-wave emission will cause them to merge in less than the age of 
the universe; and further interactions will often eject the BH/BH binaries 
from the globular, to interstel 



LT space where they merge. Simulations by 
McMillan and Portegies Zwart l3 suggest that each dense cluster will make a 
number of such BH/BH binaries, and extrapolations into the universe predict 
an event rate in WB advanced IFOs between ^ 300 and ~ 4000/yr — though 
the unccrtainti£is_are probably larger than suggested by these numbers from 
the hterature 

These event rates are very encouraging. They make it seem quite likely 
that the advanced IFOs will observe tens to thousands of BH and NS inspirals 
per year, while the initial IFOs will be lucky to observe ^ 1 per year. 

The observed inspiral waves will last for between ~ 1000 and 10, 000 cycles 
depending on the binary's masses, and will carry detailed information about 
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the binary and about general relativistic deviations from Newtonian gravity. 
This information can be extracted with good precision using the method of 
matched filters. Specifically (denoting by A^Upi Mi + M2 the binary's total 
mass and fi — M1M2/M its reduced mass): EIlO 

(i) The binary's chirp mass Mc = /i'^/^M^/^ will typically be measured, 
from the Newtonian part of the signal's upward frequency sweep, to 0.04% 
for a NS/NS binary and ~ 0.3% for a system containing at least one BH. (ii) 
//we are confident (e.g., on a statistical basis from measurements of many 
previous binaries) that the binary's spins are a few percent or less of the maxi- 
mum physically allowed, then the reduced mass fi will be measured to ~ 0.5% 
for NS/NS and NS/BH binaries, and ~ 2% for BH/BH binaries, (iii) Because 
the frequency dependences of the (relativistic) /i effects and spin effects are 
not sufficiently different to give a clean separation between n and the spins, 
if we have no prior knowledge of the spins, then the spin//i correlation will 
worsen the typical accuracy of by a large factor, to ~ 50% for NS/NS, 
~ 90% for NS/BH, and a factor ^ 1.5 for BH/BH. These worsened accuracies 
should be improved significantly (though we do not yet know hnw miir|-i) by 
waveform modulations due to spin-induced precession of the orbit EjOO, and 
even without modulational information, a certain combination of and the 
spins will be determined to a few per cent, (iv) The distance to the binary 
("luminosity distance" at cosmological distances) can be inferred, from the 
observed waveforms, to a precision ^ 3/p ^ 30%, where p = S/N is the am- 
plitude signal-to-noise ratio in the total LIGO network (which must exceed 
about 8 in order that the false alarm rate be less than the threshold for detec- 
tion), (v) With the aid of VIRGO, the location of the binary on the sky can 
be inferred, tw time of flight between the detector sites, to a precision of order 
five degrees E3. (This angular resolution is rather worse than one might ex- 
pect from a simple estimate based on time-of-flight alone. The reason is that 
because two LIGO detectors are nearly parallel, the addition of VIRGO to the 
network still leaves a near-degeneracy between polarization and time-of-flight 
information. The addition of a fourth detector with comparable sensitivity in 
Australia would decrease the error box on the sky Afi by a factor ~ 4 E3.) 

Advanced LIGO will likely produce a catalog of hundreds or thousands of 
binary inspirals and their inferred parameters; this catalog will be a valuable 
data base for observational astronomy and cosmology. 

Important examples of the general relativistic effects that can be detected 
and measured with precision, in the inspiral waves, are these: (i) As the waves 
emerge from the binary, some of them get backscattered one or more times off 
the binary's spacetime curvature, producing wave tails. These tails act back 
on the binary, modifying its radiation reaction force and thence its inspiral rate 



cutlerThorne: submitted to World Scientific on February 5, 2008 11 



in a measurable way. e3 (ii) If the orbital plane is inclined to one or both of the 
binary's spins, then the spins drag inertial frames in the binary's vicinity (the 
"Lense-Thirring effect"), this frame drageiiig causes the orbit to precess, and 
the precession modulates the waveforms c3. This precession and papdulation 
should be very strong in a significant fraction of NS/BH binaries cB. 



2.4 Tidal Disruption of a NS by a BH: Measuring the Nuclear Equation 
0/ State 

As a NS/BH binary spirals inward, its NS experiences ever increasing tidal 
forces from the BH's gravitational field (its spacetime curvature). In many 
cases these tidal forces may tear the NS apart before it begins its final, quick 
plunge into the hole's horizon. The gravitational waves from this tidal disrup- 
tion and from the termination of inspiral should carry detailed information 
about the NS's equation of state (the equation of state of bulk nuclear matter 
at ~ 10 times the density of an atomic nucleus). The disruption waves lie 
largely in the frequency band ~ 300Hz < / ^ 1000 Hz, where the wide-band, 
advanced IFOs have good sensitivity, and where narrow-band or some other 
optimized IFO can do even better. This suggests that the advanced IFOs may 
be able to extract new information about the nucleat equation of state from 
the tidal-disruption waves. A first, crude estimate suggests, for example, 
that for NS/BH binaries at 140 Mpc distance (where the event rate could be a 
few per year; see the table in Box 1), tidal-disruption observations may enable 
the NS radius R to be measured to a precision ^ 15%, by contrast with its 
present uncertainty (for fixed NS mass) of about a factor 2. From the mea- 
sured radii would follow the desired equation-of-state information. Detailed 
numerical-relativity simulations will be required to firm up this estimate, and 
will be essential as a foundation for interpreting any tidal-disruption waves 
that are observed. 

The merger waves from NS/NS binaries, by contrast with NS/BH, are 
likely to lie outside the band of good advanced-IFO sensitivity — at frequen- 
cies / ^ 1500 Hz. However the onset of NS/NS merger, triggered by a plunge 
of the two NS's toward each other, may produce a strong "cliff" in the waves' 
spectrum, in a range / ~ 400 — 1000 Hz of good sensitivity, and by measur- 
ing the cliff frequency we may strongly constrain the NS radii (assuming the 
NS masses have been determined accurately from the earlier portiop_of the 
waveform) , and hence also learn about the nuclear equation of state 



cutlerThorne: submitted to World Scientific on February 5, 2008 12 



2.5 BH/BH Mergers and Ringdown: Observing the Nonlinear Dynamics 
of Spacetime Curvature 



A BH/BH binary evolves through three epochs: a gradual, adiabatic inspiral 
driven by gravitational radiation reaction (discussed in Sec. 2.3 above); a 
merger epoch, during which the holes' spacetime geometry undergoes violent, 
highly nonlinear oscillations; and a ringdown epoch, during which the merged 
black hole undergoes damped pulsations as it settles down into its final, Kerr- 
metric, quiescent state. The inspiral waves are rather well understood (Sec. 
2.3), as are the ringdown waves; but the merger dynamics and waves are 
hardly understood at all. 

It is likely that in many cases the holes are spinning rapidly and the 
system's three angular momentum vectors (two spins and one orbital) are 
substantially misaligned. In such cases, the merger dynamics might be quite 
rich, with spin-spin and spin-orbit coupling perhaps strong enough to produce 
spin flips and complex contortions and flows of the spacetime geometry as the 
two holes convert themselv£aiato one. 

Flanagan and Hughes oEJ have made plausible estimates about the un- 
known aspects of the waves from each epoch — inspiral, merger and ring- 
down — and have used these estimates to evaluate the signal-to-noise ratios 
in LIGO and LISA for equal-mass, fast-spinning BH binaries. Figures^ and || 
show their results for the initial (LIGO-I) and advanced (LIGO-II) IFOs, and 
Fig. I (Sec. O below) shows them for LISA. [The LIGO-II results in Fig. | 
are updates 13 of the original Flanagan-Hughes results, taking account of the 
currently planned noise curve for LIGO-II, Fig. |l].] 

Figures | and | suggest that, for M ^ 24M0 (two 12Mo holes), the 
inspiral waves produce most of the signal in LIGO; for 24Mq <^ M < 2OOM0 
the merger waves dominate; and for 2OOM0 <J M <J lOOOM©, the merger 
and ringdown waves have roughly equal strength. (We say "suggest" because 
the plausible assumptions made by Flanagan and Hughes are far from secure; 
numerical-relativity simulations and/or observations are needed to firm them 
up.) For M = 20Mq (two lOAfg holes), LIGO's initial and advanced IFOs can 
detect the waves out to lOOMpc and z = 0.4 respectively, at which distance 
population synthesis calculations suggest event rates <^ 1/year and ~ 30 to 
~ 4000/year respectively (Box 1 and Sec. 2.3 above). Black-hole binaries 
more massive than M ~ 50Mq (two 25Mq holes) are unlikely to form from 
main-sequence progenitors, but could well form in dense star clusters. The 
event rates at the distances seen by the initial IFOs (^ 200 to 800 Mpc for 
masses up to ^ IOOOMq) and by the advanced IFOs (redshifts z ~ 1 to 2 
for masses up to ~ 2OOOM0) are exceedingly uncertain; but for the advanced 
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Figure 2. The inspiral, merger, and ringdown waves from equal-mass black-hole binaries as 
observed by LIGO's initial interferometers: The distance to which the waves are detectable 
(right axis) and the signal-to-noise ratio for a binary at IGpc (left axis), as funcliiops of the 
binary's total mass (bottom axis). (Figure adapted from Flanagan and Hughes e3.) 




Figure 3. (a) The signal-to-noise ratio for the inspiral, merger and ringdown waves from 
fast-spinning, equal-mass black-hole binaries at a fiducial distance of IGpc, as observed 
by advanced (LIGO-II) interferometers, (b) The cosmological redshift out to which these 
waves can be seen, assuming a cosmology with Hubble expansion rate Ho = 65 km/s/Mpc, 
cold-matter density 0.4 of that required to close the universe Qm = OA, andivacuum energy 
density (cosmological constant) 0.6 of closure, Qa = 0.6. (Figure from Ref. Q, basedreri the 
LIGO-H wide- band noise curve in Fig. hi and the analysis of Flanagan and Hughes Ej ) 
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IFOs they could be large. 



2.6 Supernovae and Accretion-Induced Collapse of White Dwarfs 

Type-II supernovae are triggered by the violent collapse of a stellar core to 
form a NS or BH. Despite decades of extensive modeling, we still understand 
only poorly the details of the collapse and how it triggers the explosion, and 
the birth throes of the NS or BH. The optical display brings little information, 
as the stellar material is so optically thick that the display is produced hours 
after the collapse and from radii ~ 10,000 times larger than the NS or BH. 
Neutrinos bring information about the core's temperature during the first few 
seconds after the collapse, and gravitational waves could reveal the nonspher- 
ical dynamics of the core's mass distribution on millisecond timescales and 
longer, and the nature of the core's final state: BH or NS. 

Numerical models of supernovae suggest that, even if it is slowly rotat- 
ing, a newborn proto-ncutron-star will be convectively unstable, and that the 
gravitational waves from the convective overturn in the first ~one second of 
the star's life may be detectable tkroughout our galaxy and its orbiting com- 
panions, the Magellanic Clouds 113. Although the supernova rate is low in 
our galaxy and its companions (< l/30yrs), one observed event could be very 
valuable scientifically: The bulk of the supernova's neutrinos are thought to 
come from the same convecting material as produces the gravitational waves, 
so there should be correlations between the neutrinos and the waves, which 
could teach us much about the proto-neutron-star's dynamics. 

Observations show that many type-II supernovae produce neutron stars 
and give them kicks of magnitude as large as ~ 1000 km/s. Indeed—about 
half of all radio pulsars are born with a kick larger than 500 km/s cj. The 
observed kicks suggest that at least some newborn proto-neutron-stars may 
be strongly asymmetric, perhaps due to fast rotation, and therefore could 
produce significant gravitational radiation. Rapidly rotating proto-neutron- 
stars may also be produced by the accretion-induced collapse of white dwarf 
stars (AIC) in some cases, depending on the white-dwarf composition, central 
density and accretion rate E3. 

If the newborn proto-neutron-star in an AIC or supernova is spinning 
so fast that it hangs up centrifugally at a radius large compared to that of 
the final NS, then it may be dynamically (or at least secularly) unstable to 
deforming into a bar-shaped object that tumbles end-osier-end, emitting grav- 
itational waves in LIGO's band of good sensitivity l3. Recent simulations 
suggest that the bar will be long-lived, rather thari just wrapping itself up 
into an axisymmetric shape and disappearing Cj'eZI. The simulations show 
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that the waves from such a bar may sweep upward in frequency, due to a 
gradual shrinkage of the proto-neutron-star, or may S|S4£eep downward due 
to development of "Dedekind-like" internal circulation e3. The discovery of 
the waves from such a proto-neutron-star and observations of their frequency 
evolutions (which should mirror the bar's evolution) would teach us much. 
Though the strengths of the waves and the best signal processing techniques 
and thresholds are all ill-understood, a rough estimate lij suggests that the 
advanced IFOs' range for detection might be the distance of the VIRGO clus- 
ter of galaxies (about 15 Mpc) and conceivably significantly larger, suggesting 
event rates that could be some per year but might be far less. 

If the advanced WB IFOs detect no gravitational waves from a supernova 
at distance r, during a time T preceeding the beginning of the optical outburst 
or over a time T during a neutrino outburst, then one can thereby place a limit 
on the emitted gravitational-wave energy. The limit (as devised by FinnEJ) 
is based on the statistics of a cross correlation of the outputs of the two 4km 
LIGO interferometers, over the chosen time T and over some frequency band 
near the minimum of the advanced LIGO noise curve. In the band from ~ 100 
to 300 Hz this limit is AE^gw <, O.O5(r/15Mpc)2(T/lh)i/2j\^0. In the case of 
a neutrino-emitting supernova in our own galaxy, with T ^ 1 sec, this limit 
is very impressive: 10~^Mq. 



2. 7 Spinning Neutron Stars 

As a newborn neutron star settles down into its final state, its crust solidifies 
(crystalizes). The solid crust will assume nearly the oblate axisymmetric 
shape that centrifugal forces are trying to maintain. Thereafter the crust has 
a preferred shape (oblate about some preferred axis). One can define the 
NS's deformation ellipticity e^j as the residual ellipticity the NS would then 
maintain (due to induced crustal shear stresses) if the NS were spi].n .down to 
zero frequency without the crust breaking or otherwise relaxing e1[e3: 

.,.7xl0-(ilif (1) 

where v^ci is the spin frequency for which the crust is (most) relaxed. Again, 
this is the level of ellipticity that is "frozen into" the crust, as opposed to the 
ellipticity that arises from centrifugal forces and so "follows" the star's spin . If 
the star's angular momentum J deviates from the crust's preferred symmetry 
axis, the NS will wobble as it spins, with small "wobble angle" 0^,. The NS will 
emit GW's with h cx ^^e^ at frequency few = /rot + /prcc, where the body 
frame precession frequency fp^ec [= ed/iot(-fcrust/^)] is orders^Ljnagnitude 
smaller than /lot; thus few is just slightly higher than /rot E3H. The NS 
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may also deviate slightly from axisymmetry about its principal axis; i.e., it 
may have a slight ellipticity eg = (/i — l2)/l3 in its equatorial plane. It will 
then emit GW's with h oc Ce at twice the rotation frequency, few = 2/rot 
The gravitational ellipticity e referred to in Fig. |l| and in Sec. 2.2 is e — ee 
for deformations around the principal axis, and e = 29wed for deformatins in 
which J is inclined to the preferred symmetry axis, producing wobble. 

How large will be eg and O^ed'^ Ushomirsky et al. 113 have shown that, 
for typical NS parameters, ee < 5 x 10~^((Tma£c/10~^), where amax is the 
(highly uncertain) crustal yield strain, estimated to be somewhere in the range 
10~^ < amax < 10^^. Thus rough upper limits are ee < 10~^ and [from 
Eq. 1^ Ow£d < 10~^. These upper limits are for ellipticities that are basically 
supported by crustal shear stresses, but there are other possible sources of 
ellipticity. For example, inside the superconducting cores of many neutron 
stars there may be trapped magnetic fields with mean strength Score ^ 10^'^G 
or even lO^^G. Because such a field is actually concentrated in flux tubes 
with B = Bait ~ 6 X lO^^G surrounded by field-free superconductor, its 
mean pressure is pB = i?coroScrit/87r. This pressure could produce ellipticity 
Ec Ed ~ Pb/p ^ lO^^Bcoro/lO^^G (whcrc p is the core's material pressure). 
We are extremely ignorant, and correspondingly there is much to be learned 
from searches for gravitational waves from spinning neutron stars. 

From this discussion we see that the range of plausible gravitational ellip- 
ticities for a spinning neutron star is e < 10"^, though they may get as large 
as 10~^. The strengths of the waves produced by such ell ipti cities are shown 
in Fig. |l| and the waves' detectability is discussed in Sec. 2^. 



Neutron stars gradually spin down, due in part to gravitational- wave emis- 
sion but perhaps more strongly due to electromagnetic torques associated with 
their spinning magnetic fields and pulsar emission. This spin-down reduces 
the strength of centrifugal forces, and thereby causes the star's oblateness to 
decrease, with an accompanying breakage and resolidification of its crust's 
crystal structure (a "starquake") M. In each starquake, 0^, ee, and e^ will 
all change suddenly, thereby changing the amplitudes and frequencies of the 
star's two gravitational "spectral lines" / — 2frot and / = /rot + /prcc- Af- 
ter each quake, there should be a healing period in which the star's fluid 
core and solid crust, now rotating at different speeds, gradually regain syn- 
chronism. By monitoring the amplitudes, frequencies, and phases of the two 
gravitational- wave spectral lines, and by comparing with timing of the elec- 
tromagnetic pulsar emission, one might learn much about the physics of the 
neutron-star interior. 

A spinning NS can radiate not only due to built-in asymmetries, as dis- 
cussed above, but also due to a gravitational-radiation-reaction-driven insta- 
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bility first discovered by Oiandrasekhar cB and elucidated in greater detail 
by Friedman and Schutz E3. In this "CFS instability" , density waves travel 
around the star in the opposite direction to its rotation, but are dragged for- 
ward by the rotation. These density waves produce gravitational waves that 
carry positive energy as seen by observers far from the star, but negative en- 
ergy from the star's viewpoint; and because the star thinks it is losing negative 
energy, its density waves get amplified. This intriguing mechanism is similar to 
that by which spiral density waves are produced in galaxies ^Allhough the CFS 
instability was once thought ubiquitous for spinning stars Ej'lZI, we now know 
that neutron-star viscosity will kill it, stabilizing the star and turning off the 
waves, when the stan's temperature is above some limit ~ lO^^K L3 and below 
some limit ^ lO^KEj; and correspondingly, the instability should operate only 
during the first few years of a neutron star's life, when lO^K <,T <^ 10^°K, and 
only for NS spins, greater than ^ 90% of the upper limit set by mass-shedding 
at the equator E2l. 

Shortly after GRlSjthere was a surge of renewed interest in the CFS insta- 
bility when Andersson Eil and Friedman and Morsink o showed that r-modes 
(driven by Coriolis forces) are unstable al_all spin values (again, in_the absence 
of viscosity), and then Lindblom et al. S and Andersson et al. l3 estimated 
that in a newborn NS, the instability would not be damped by shear or bulk 
viscosity, but instead would rapidly spin down the NS to v < 100 Hz. How- 
ever, since then, a host of other damping-jpechanisms have been suggested, 
includinsp. (i) wind-up of magnetic fields E^, (ii) dissipation at the crust-core 
bondaryEa, (iii) bulk viscous damping from weak interactions involving hyper- 
ons in the inner core, which is ppedirted to kill the instability unless the star is 
less masive than about 1.2Mq Ellea; and (iv) nonlinear mode- mode coupling, 
specifically a 3-mode interaction, which is predicted to damp out r-modes-at a 
maximum dimensionless amplitude ~ 10~'^(spin frequency/lOOOHz)'^/^ e3. It 
therefore appears that r-modes will not be significant GW sources, after all. 



2.8 Gamma Ray Bursts 

Cosmic gamma-ray bursts, observed ~ once per day by detectors on space- 
craft, are thought to be produced by extremely hot shocks in a relativistic 
fireball, which is ignited by rapid inflow of gas onto a newborn black hole l3. 
The fireball is so dense that electromagnetic observations cannot probe the 
BH's vicinity; the observed gamma rays are produced at radii ^ lO^'^ to 10^^ 
cm — millions to billions of times larger than the hole itself. The only hope for 
probing the fireball's trigger is via gravitational waves: The BH and its accre- 
tion flow are thought to form violently, by the collapse of a massive star (per- 
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haps initiated by merger with a companion) ["a hypernova"], or by the merger 
of a binary made of compact objects: ajNS/NS binary, a BH/NS binary, a 
BH/white-dwarf or a BH/He-core binary E2l. Each of these gamma-burst trig- 
gers should emit strong gravitational waves that carry detailed information 
about its source. 

The distribution of 7-ray bursts is bimodal: there are long bursts (du- 
ration ^ 2 s) and short bursts (<; 2 s). Beppo-Sax has detected afterglows 
from several long bursts, allowing their redshifts to be measured. No red- 
shifts are known for stuirt bursts, but their distance can be estimated from 
the (V/Vmax) statistic lid (assuming a standard profile, and that their event 
rate tracks the star formation rate). The best event rate estimates (by M. 
Schmidt lid) currently are: 1/yr to D = 17O(0o/O.l)^^^ Mpc for long bursts, 
and 1/yr to D = 25O(0o/O-l)^^^ Mpc for short bursts. Here the gamma-rays 
are assumed beamed into solid angle 2tt9q. Clearly, these rates (for 9o = 0.1) 
are quite consistent with the estimated NS-NS merger rates given in Box 1, 
and marginally consistent with the BH-NS rate. 

If BH-NS mergers are indeed responsible for either long or short bursts, 
the (WB) LIGO II detection rate would be ~ 15 — 50/yr. Of course, these 
closest sources will generallly not be seen in 7-rays, since they are typically 
beamed away from us; the closest burst for which we see the 7-rays (in 1 yr) 
is at I? = 1.0(1.5) Gpc for long (short) bursts. Could these be detected in 
coincidence with GW inspiral events? Box 1 states that LIGO II (WB) could 
detect BH-NS (NS-NS) inspirals to 0.65 (0.30) Gpc. But those distances 
assume a random orientation of the source, and assume that S/N > 8 is 
required for confident detection. If the 7-rays are beamed perpendicular to 
the orbital plane, then in cases where we detect the burst, we are also at 
optimal "viewing angle" for the GW's. The amplitude h of the GW's beamed 

"up the axis" are ~ 1.58 larger than the rms value (averaged over the 

sky, as seen from the source). Also, the near-coincident detection of a 7-ray 
burst in time (within several seconds of each other) and on a consistent area of 
the sky, effectively reduces by pother factor 1.5 the GW S/N required for 
detection with 99% confidence E3. Thus the distances out to which LIGO I and 
II could detect binary mergers associated with 7-ray bursts are ~ 2.4 times 
farther than the distances given in Box 1; for LIGO II: 720 Mpc for NS/NS 
and 1.55 Gpc for NS/BH mergers. Thus, if NS/BH mergers are responsible for 
long and/or short bursts, with LIGO II one might expect ~ 2/yr simultaneous 
detections in both GW's and 7-rays. If NS/NS mergers are responsible, one 
might expect to wait ~ 5 yr. 

M. van PuttenO points out that in the above 7-ray burst scenarios, the ac- 
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cretion torus itself might develop large non-axisymmetries and so be a strong 
source of GW's. He develops a picture of long bursts as due to "suspended" 
accretion around a rapidly rotating BH, with most of the BH's angular mo- 
mentum getting transferred to the torus via Maxwell stresses (similar to the 
Blandford-Znajek effect). If the torus develops "lumps" , van Putten estimates 
that ~ 10% of the BH's spin energy could get radiated away in GW's-so 
Egw ~ Q-ZMqIMbh I^^Mq). He predicts a sinusoidal waveform with fre- 
quency / ~ 1-2 kHz that wanders slowly in time — which is simple enough 
to allow near-optimal data analysis. Would this be detectable? The rms S/N 
would then be (taking z \ and approximating / and Sh{f) as constants 
over the signal band, for simplicity), 

S _ \/2Egw 

yielding 6*//^ = 1.4 for Egw = 0.3Mq, D = 150 Mpc, / = lO^Hz, and 
hni= V^fShif)) = 8 X 10-22 (LIGO II WB value at 1 kHz). This does 
not seem especially promising, but given the uncertainty in Eqw and D (for 
the closest source in one year), it seems worth looking for. Moreover, one 
could improve the sensitivity in LIGO-II by a factor of a few by reshaping 
the noise curve so it is lower at / '--^ 1 kHz, at the price of higher noise at 
lower frequencies (cf. the difference between the NB and WB LIGO-II curves 
inFig.|l|). 

If gravitational waves are detected from one or more gamma-burst trig- 
gers, the waves will almost certainly reveal the physical nature of the trig- 
ger. Moreover, by comparing the arrival times of the gravitational waves 
and the earliest gamma rays, it should be possible to measure the rel- 
ative propagation speeds of light and gravitational waves to an accuracy 
~ 1 sec/lO^" yr- lO^^^. 

If no gravitational waves are detected from any individual gamma burst, 
then (as Mohanty, Finn and Romano lU have shown) the correlation between 
gamma bursts and gravitational waves might nevertheless be established by 
statistical studies of the advanced IPOs' gravitational-wave data in narrow 
time windows preceeding the gamma bursts. 

2. 9 Sources of Stochastic Background 

Searches for a GW stochastic background with Parth-based detectors begin 
with the assumptions that 1) the GW background will certainly be smaller 
than the instrumental noise at all frequencies, and 2) the instrumental noise 
level will not be known sufhciently well a priori that one can simply search 



cutlerThorne: submitted to World Scientific on February 5, 2008 20 



for "excess noise" in each instrument. Thus, stochastic background searches 
with Earth-based detectors are based on cross-correlating the outputs of two 
highly independent detectors. (The two 4 km LIGO detectors at Hanford, 
Washington and Livingston, Louisiana give the best sensitivity). 

The most plausible sources of a stochastic gravitational- wave background 
in LIGO's frequency band are processes in the very early universe. The 
current best limit on the strength of such waves is 6 x 10^^, where 

= J il{f)d{\n f) is the integrated GW energy density from all frequency 
bands. A wave energy larger than this would have caused the universe 
to expand too rapidly through the era of primordial nucleosynthesis (uni- 
verse age ^ a few minutes), thereby distorting the universal abundances 
of light elements away from their observed values. Other known bounds 
on fl{f) are: 1) /i^QQf2(/ ^ lO^^Hz) < 10"^ from pulsar timing; and 
2)hl„on{f) < 7 X 10-"(//^o)"' for 3 x IQ-^^hiooH^ <f< xlQ-^^hwo^z, 
from microwave background anisotropy measurements. Here ft-ioo is the Hub- 
ble constant in units of 100 km/s/Mpc. 

The LIGO H interferometers would improve the current limit at / ~ 40 
Hz by a factor ^ 10^, to ^l{f ~ 40Hz) ~ 5 x 10~^ — an improvement enabling 
LIGO to test a number of current speculations about the very early universe. 
Inflationary models of the early universe predict that vacuum fluctuations, 
created in the Planck era when the universe was being born, should have been 
parametrically amplified, during the first ^ 10~^^ sec of the universe's life, to 
produce a stochastic gravitational wave background in the LIGO band. Un- 
fortunately, if "standard" inflation theory is correct, tiifin the amplified waves 
are much too weak for LIGO to detect, fl ^ 10^^^ L3. The most plausible 
modifications of standard inflation push downward from this E3, but spme 
less plausible modifications push it upward, to the point of detectability 119. 

The first, tentative efforts to combine superstring theory with inflationary 
ideas have produced a new description of the very early universe called the 
"pre-big-bang model" , in which string effects cause the gravitational- wave 
spectrum to rise steeply at high frequencies — most likely jat frequencies 
above LIGO's band, but quite possibly in or below that bandllH. The result 
could be waves strong enough for LIGO to detect. A non-detection would 
significantly constrain the pre-big-bang model. 

There are a wide variety of postulated mechanisms that could have pro- 
duced strong gravitational waves, with wavelengths of order the horizon size, 
at various epochs in the very early universe. Those waves (if any) produced at 
(universe age)~ 10~^^ sec, corresponding to (universe temperature T) ~ 10^ 
GeV, would have been redshifted into the LIGO band today and might be 
detectable. The temperature (energy) region 10^ GeV is tera incognita; 
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LIGO's advanced detectors will provide our first opportunity for a serious 
experimental exploration of it. Among the speculated wave-generating mech- 
anisms that could operate there, and that LIGO could constrain (or discover!), 
are these: 

• A first-order phase transition in the states of quantum fields at T '--^ 10^ 
K. Such a phase transition would nucleate bubbles of the new phase that 
expand at near the speed of light and collide to produce gravitational 
waves; and their collisions would also generate turbulence that radiates 
waves. If the transition is strongly first order—the waves would be strong 
enough for LIGO's advanced IFOs to detect 113. 

• Goldstone modes (coherent, classical excitations) of scalar fields that arise 
in supersymmetric and string theories. If strongly excited, these modes 
will entail coherent Hqpis of energy that radiate gravitational waves strong 



• Coherent excitations of our 3+1 dimensional universe, regarded as a 
"brane" (defect surface) in a higher dimensional universe. The excita- 
tions could be of a "radion" field that controls the size or curvature of 
the additional dimensions, or they could be of the location and shape of 
our universe's brane in the higher dimensions; in either case, if there is 
an equipartition of energy between these excitations, in the very early 
universe, and other forms of energy, then the excitations will produce 
gravitational waves easily strong enough for detection by LIGO's ad- 
vanced IFOs. LIGO would thereby probe one or two additional dimen- 
sions of size or curvature length ~ 10"^" ~ 10"^'^ mm; by contrast, LISA's 
lower-frequency observations would probe lengths ~ 1 - 10~^ mm. If the 
number of extxajdimensions is larger than 2, the probes reach to much 



Cosmic strings (not to be confused with superstrings) , produced in the 
early universe, were once regarded as candidates for seeding galaxy formation, 
but recent cosmological observations have ruled them out as seeds. Neverthe- 
less, it remains possible that a network of vibrating cosmic strings too weak to 
seed galaxy formation was formed in the early universe. LIGO can search for 
the presence of such a network in two ways: (i) Via the stochastic background 
of gravitational waves that the strings' vibrations produce; this background 
would be strong enough fjor the advanced IFOs to detect if the strings' mass 
per unit length is ^ 10~^. e3 (ii) Via occasional, individually detectable, strong 
non-Gaussian bursts ( Jlspikes" ) of gravitational waves produced by kinks or 





cusps along the string 
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Cusps are points on the string that, at one instant, move at the speed of 
hght; kinks typically arise when two strings cross each other and interconnect. 
We now briefly discuss the associated GW spikes as analyzed by Damour 
and Vilenkin EJ. These signals have no free parameters, except the overall 
amplitude and the arrival time tc'- h{t) oc |i — icP^'^ for cusps and h{t) oc 
\t — tc\^/'^ for kinks. The time spent in the detector band is ~ l//c where 
/c is the center of the band. Matched filter searches should therefore be 
trivial, with coincidence between detectors (and other vetoes) used to remove 
instrumental bursts. Cusp bursts are more highly beamed than kink bursts, 
have higher amplitude, and are more detectable (the higher amplitude more 
than compensates for the decreased probability that any given cusp burst is 
beamed in our direction). Cusp bursts should be detectable for a large range 
of string tension /i not otherwise constrained: LIGO II should detect cusp 
bursts for /x ^ 10~^^, while LISA could detect n ^ 10"^'^. Burst searches 
thus provide a stronger limit on cosmic strings than searches for a stochastic, 
Gaussian background. 

Each of these cosmological speculations is plausible, though not highly 
likely. (For further explication of most of them, we highly recommend the 
review by Allen Efl.) Perhaps their greatest value is to remind us of how 
terribly ignorant we are of physics and astrophysics in the domain that LIGO's 
advanced IFOs will probe. Our ignorance may well be even greater than that 
of the pioneering radio astronomers of the 1930s and X-ray astronomers of 
the 1960s; and as there, so also here, the first waves to be discovered may 
well be from sources that were previously unknown. LIGO II and its partners 
could bring us a revolution of insights into the universe comparable to the 
revolutions wrought by radio and X-ray astronomy. 

3 Low-Frequency Sources Sought by LISA 

3.1 LISA'S Noise Curve and Conventions for Wave Strengths 

In discussing low-frequency sources, we shall place them in the jContext of the 
sensitivity of LISA, the Laser Interferometer Space Antenna E3. As is dis- 
cussed in Schutz's chapter in this volume □, LISA consists of three drag-free 
spacecraft at the corners of a 5-million-kilometer triangle. The spacecraft 
track the distances to each other using laser beams, thereby searching for 
gravitational waves in the 10~^-1 Hz low-frequency band. LISA is tentatively 
scheduled for launch in 2011 as a joint European- American mission with a 
nominal lifetime of 5 years, though it could well continue to collect GW data 
much longer than that. By the time LISA begins its GW searches, the second- 
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Figure 4. HSA's sky- and polarization-averaged rms noise h„{f) in a bandwidth equal to 
frequency, compared with the strengths of the waves from several low-frequency sources. 
The wave strengths are plotted so the height above the noise curve is the S/N for a wave 
search using optimal signal processing; see text for further detail. For frequency-sweeping 
waves (arrowed curves) the dots are, from left to right, the signal strength and frequency 
1 year, 1 month, and one day before the end. The thick-dashed curve is a background of 
waves from WD- WD binaries that are so numerous they cannot be removed from the data, 
and so constitute a noise when searching for other waves. 



generation earth-based detectors (e.g. the advanced LIGO-II IFOs) described 
in the last section could be nearing the end of their observations, and a tran- 
sition to a third generation might be near at hand. 

Figure ^ depicts LISA's noise and compares it with the wave strengths for 
various low-frequency sources, in much the same manner as Fig. ^ depicted 
LIGO's noise and compared it with high-frequency wave strengths. However, 
there are important differences in conventions between the LISA Fig. ^ and 
the LIGO Fig. |l|: 

First: For LIGO we plotted vertically the square root of the spectral 
density of the fractional arm-length difference, h{f) = ^/S'^j^JlU) (noise per 
root Hz) for a single 4 km IFO. For LISA, by contrast, we plot the rms 
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GW strain noise in a bandwidth equal to frequency A/ = /, averaged over 
the sky and over GW polarizations, for the full 3-spacecraft LISA system: 



hn — y f^h^if), where S^^{f)Js the noise spectral density inverse-averaged 

over the sky and polarizations eII. (For LIGO the analogous quantity is /i„ = 
^/5jh, where the -\/5 accounts for the average over the sky and polarizations.) 
In many discussions of LISA a different quantity is used to characterize the 
noise: the amplitude sensitivity to periodic sources for a one-year integration 
time and a signal-to-noise ratio of 5, averaged over the sky and over source 
polarizations. That quantity, /i|n5 lyri is related to our rms noise hn by 



_ V/ X lyr SA 



Second: Our wave-strength convention for LISA (Fig. is also different 
from that for LIGO (Fig. 0): The height of a source point or curve above 
the LISA noise curve is equal to the signal-to- noise ratio S/N in a search for 
the source's waves using optimal signal processing. By contrast, for LIGO 
(Fig. |l]), the height above the noise curve is the signal-to-threshold ratio S/T, 
with the threshold being that at which the false alarm probability is one per 
cent when using the best currently known, practical data analysis algorithm 
(typically T ^ 57V). The difference in convention is motivated by the fact 
that we know little as yet about practical LISA data analysis algorithms and 
thus don't know with much confidence how to set thresholds, whereas LIGO 
data analysis is fairly mature. 

For LISA, when evaluating the S/N^s that determine the height of the 
source points in Fig. ^, we assume for periodic sources (e.g. white dwarf 
binaries) an integration time of five years (the nominal mission lifetime), and 
for sources with a large frequency sweep, we assume at a given frequency / 
an integration time that is the shorter of the time until the signal stops, and 
the time to sweep through a frequency band A/ = /. 

Notice that LISA's minimum noise /i„ ~ 10~^^ is nearly the same as 
that of LIGO's first interferometers (Fig. with /i„ = V^fh), but at 100,000 
times lower frequency /. Since the waves' energy flux scales as f^h^, this 
corresponds to 10^° better energy sensitivity than LIGO-I. 

LISA can detect and study, simultaneously, a wide variety of different 
sources scattered over all directions on the sky. The key to distinguishing 
the different sources is the different time evolution of their waveforms. The 
key to determining each source's direction, and confirming that it is real and 
not just noise, is the manner in which its waves' amplitude and frequency 
are modulated by LISA's complicated orbital motion — a motion in which the 
interferometer triangle rotates around its center and around the normal to the 
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ecliptic plafle once per year as the spacecraft orbit the sun; see the discussion 
by Schutz □ elsewhere in this volume. Most sources will be observed for of 
order a year or longer, thereby making full use of these modulations. 

Our understanding of sources and data analysis for LISA has evolved con- 
siderably since GR15. The principal sources are (i) short-period stellar-mass 
binaries — both galactic and extragalactic, (ii) the inspiral and merger of 
massive (~ IO^Mq) BH binaries, (iii) the inspiral and capture of stellar-mass 
compact objects (white dwarfs, NS's and BH's) into massive BH's, (iv) the 
gravitational coUape of supermassive stars, and (v) a stochastic GW back- 
ground from the early universe. We discuss these in turn in the following 
sections. 

3.2 Short- Period Binaries 

Unlike the case for ground-based detectors, there are known GW sources that 
LISA is certain to detect at its planned sensitivity: short-period binary stars 
in our galaxy. There are currently about a dozen galactic binaries with GW 
frequencies above 0.1 paHz that LISA would be able to detect and study in 
detail. These include: E3 (i) White-Dwarf / white-dwarf (WD/WD) binaries; 
the most detectable known example is WD 0957-666 (shown in Fig. ^), with 
masses 0.37 and O.32M0, distance from Earth 100 pc, and time to merger 
2 X 10* years, (ii) Am CVn stars — a white dwarf (WD) that accretes from 
a low-mass helium-star companion; the most detectable known example is 
RXJ1914-I-245 (shown in Fig. |), a 0.071/© hehum star orbiting a 0.6Mq 
WD at lOOpc distance from Earth, (iii) Low-mass X-ray binaries — a NS that 
accretes from a low-mass companion; the most detectable known example is 
4U1820-30 (shown in Fig. |), a < O.IA/© star orbiting a IAMq NS at 8100 
pc from Earth. These sources will provide an important test that LISA is 
functioning as expected (as well as a test of the general relativistic theory of 
GW emission, in the weak-gravity limit). 

These optically-known binaries are just the tip of the iceberg; there are 
an estimated 10®~^ galactic binaries with GW frequencies / > 0.1 mHz (most 
are WD/ WD). At frequencies below ^ 2 mHz, there are many binaries per 
resolvable frequency bin A/ — 1/Tobs ^ lO^^Hz — too many to be fit and re- 
moved from LISA'S data, so they constitute a source of "confusion noise" that 
will dominate over the instrumental noise at these frequencies. Thus LISA 
truly suffers from an "embarrassment of riches." But there is also some science 
to be extracted from this confusion noise: because the WD/WD background 
is anisotropic, the confusion noise level is modulated by LISA's time-varying 
antenna pattern on the sky. From this one can extract with few-percent ac- 
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curacy the m = I and m — 2 moments_(in spherical coordinates centered on 
the Sun) of the WD/WD distribution El. 

Above ~ 2 mHz the density of sources thins to less than ~ 0.5 per fre- 
quency bin, and one expects to be able to measure individually all these 
binaries, about 10* of them, and remove them from the data. The confusion 
noise therefore drops dramatically above ^ 2 mHz — to a level set by the 
extragalactic binary background (about 10 times smaller than the galactic 
background) — and crosses below the instrumental noise curve; see Fig. ^ 

Although it should be possible to remove the binaries above 2 mHz from 
LISA'S data so as to search there for other waves (e.g. from inspiral of small 
mass BH's into supermassive BH's; cf. Fig. U), no effort has yet been put 
into the task of developing practical data analysis algorithms to perform the 
removal. Designing and optimizing such algorithms is an important challenge 
for the next few years. 

For the 10^ highest-frequency binaries, one should be able to measure 
the binary's rate of change of frequency / (due to inspiral driven by GW emis- 
sion). From the measured / and / one can infer the binary's chirp mass Mc, 
and from Mc and the measured amplitudes of the binary's two polarizations 
once can infer its inclination to the line of sight and its distance. For binaries 
with / ~ 10~^ Hz, one ejspects to measure the angular location to '--^ 1° lj 
and the distance to ^ 2% til. Thereby one can obtain a 3-D map of the galaxy 
with distances unsullied by problems of dust, etc. 



3.3 Mergers of Supermassive Black Holes 

The merger of two supermassive black holes (SMBHs; lO'' — lO^Mo) at z = 1 
could be detected by LISA with S/N of thousands; see Fig. |. This S/N level 
is at least 10 times greater than we can expect with ground-based detectors 
(Figs. ^ and ||) — so high that the inspiral portion of the waveform could 
typically yield the binary parameters to the following accuracies: individual 
masses to AM^/M^ ~ 10^*, the so-called spin-orbit parameter /3 to Af3 ~ 
10~^, the distance to ADUJ ~ 10"'^ — 10~^, and the position of the source 
on the sky to ~ 10~* sr l£3. Knowing the binary's initial conditions with 
such high accuracy, it should be possible to make high-precision comparisons 
of numerical relativity simulations with the observed merger waveforms — 
thereby extending LIGO's exploration of the nonlinear dynamics of curved 



spacetime (Sec. 2.5) into a domain of far higher accuracy. For example, as 
the gravitational waves from a binary's merger depart from their source, the 
waves' energy should create (via the nonlinearity of EinsteiiiIsJB,eld equations) 



a secondary wave called the "Christodoulou memory" EIHE, which LISA 
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Figure 5. The waves from equal-mass, supermassive black-hole binaries as observed by LISA 
in one year of integration time. Plotted horizontally is the binary's total mass M multiplied 
by 1 plus the binary's cosmological redshift z. The solid, dashed, and tight-dotted curves 
refer to the left axis, which shows the estimated S/N for LISA'S observations when the binary 
is 1 Gpc from earth (redshift z ~ 0.2). The wide-spaced dots are curves of constant binary 
mass M, for use with the right axis, which shows the luminosity distance and redshift to 
which the binary can be detected, with S/N = 5 (look at the intersection of a wide-dotted 
curve with a solid, dashed, or tight-dotted curve and then project horizontally onto the 
right axis). (The assumed cosmology is Ho = 65 km/s/Mpc, = 0.6, Qm = 0.4.) The 
bottom-most curves are the signal strengths after one year of signal integration, for BH/BH 
binaries[4p years and 100 years before their merger. (Figure adapted from Flanagan and 
Hughes c3.) 



should easily be able to measure, while LIGO might not quite see it. 

Unfortunately, it is far from obvious whether the event rate for such 
SMBH mergers will be interestingly high: By contrast with advanced LIGO 
IFOs, where Kalogera and others estimate tens to thousands of BH/BH merg- 
ers per year (Box 1), Haehnelt i3 and others estimate for LISA's SMBH 
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mergers out to z = 1 a rate in the range 10""'^ — 10^ per year. Even if the 
coalescence rate is only 0.1/yr, then LISA should still see ~ 3 BH/BH bina- 
ries with 3OOOM0 <^ M <^ IO^Mq that are ~ 30 years away from their final 
merger. These slowly inspiraling binaries should be visible, with one year of 
integration, out to a redshift z ~ 1 (bottom part of Fig. ^). 

Because the SMBH- merger S/N is so high at z = 1, LISA could detect 
such mergers to redshifts z ^ 1. SMBH's are not expected to form in the 
universe at redshifts any larger than z ~ 30, and perhaps somewhat less than 
this. Figure || shows that the merger waves can be detected to distances well 
beyond z ~ 30 for masses in the range M ~ 10'' to ^ IO^Mq. (Below a few 
times 1O^M0 a crucial role in making this possible is played by the redshift 
itself, which lowers the waves' frequency, bringing them into LISA's band of 
good sensitivity; see Fig. ^.) 

In the standard picture of bottom-up structure formation, large galax- 
ies are formed by the merger of protogalaxies (which are themselves formed 
from the mergers of smaller protogalaxies). If the the merging protogalaxies 
contain IO^Mq BH's at their centers, dynamical friction will bring the BH's 
to the common center, plausibly (though not assuredly) to the point where 
gravitational radiation takes over and drives the BH's to merger in less than 
a Hubble time. Based on this picture, Phinney estimates that the merger 
rate for IO^Mq BHjsjvould be ~ 150Q/yr at z ~ 7, ~ 50Q/yr at z ~ 10, and 
~ 1/yr at z '--^ 20 E3lZI — rates comparable to those of advanced LIGO IFOs 
(Box 1). Thus LISA has the possibility to reach back and shed light on the 
earliest epochs of galaxy and SMBH formation 



3.4 Inspirals of ^ 1 ~ lOOOM© Compact Objects into Massive BH's 

LISA has the potential to observe the final inspiral waves from compact ob- 
jects (WD's, NS's and BH's) of any mass fj, ^ IMq spiraling into central 
bodies of mass 3 x lO^M© ^ M ^ 3 x lO'^M© out to 3Gpc. Figure | shows 
the example of a IOMq black hole spiraling into a IO^Mq SMBH at 1 Gpc dis- 
tance. The inspiral orbit and waves are strongly influenced by the SMBH spin. 
Two cases are shown l3: an inspiraling circular orbit around a non-spinning 
hole, and a prograde, circular, equatorial orbit around a near maximally spin- 
ning hole {a/M = 0.999 in Kerr-metric language). For each case three dots 
are shown on the signal curve — from left to right the signal strength and 
frequency one year, one month, and one day before the small BH plunges into 
the horizon of the SMBH. Table |l| shows, at each of these times, the ratio 
of the orbital circumference to the SMBH horizon circumference, C/Ch, the 
number of GW cycles left until plunge, A'cyc, and the signal to noise in a 
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Table 1. The waves from a IOMq in circular equatorial orbit around a 10 SMBH. 



SMBH Spin 


Time to Plunge 


C/Ch 


^ *cyc 


S/N 


0.999 


1 yr 


6.80 


185,000 


200 


0.999 


1 mo 


3.05 


40,900 


60 


0.999 


1 day 


1.30 


2,320 


5 





1 yr 


9.46 


85,000 


60 





1 mo 


7.14 


9,650 


70 





1 day 


6.22 


366 


20 



bandwidth equal to frequency (or, if that band extends beyond plunge, then 
the total signal to noise in the remainder of the waves), S/N. 

The numbers in Table |l| are remarkable: The small BH spends its last year 
in the immediate vicinity of the SMBH horizon, emitting ~ 100, 000 cycles 
of waves. These waves should_carry, encoded in themselves, high-accuracy 
informatioLD^about the SMBH l3: 

Ryan C3 has shown that, when any small-mass, compact object spirals 
into a much more massive, compact central body, the inspiral waves will 
carry a "map" of the massive body's external spacetime geometry. Since 
that geometry is uniquely characterized by the values of the body's multiple 
moments, we can say equivalently that the inspiral waves carry, encoded in 
themselves, the values of all the body's multipole moments. By measuring 
the inspiral waveforms and extracting their map (i.e., measuring the lowest 
few multipole moments), we can determine whether the massive-central body 
is a black hole or some other kind of exotic compact object e3. The huge 
number of wave cycles, in the example in Table |^, implies the possibility of a 
high-precision map. 

While the S/N values in Table |l| are large enough (up to 200) to seem 
encouraging for detecting these waves and extracting the map that they carry, 
there are serious worries. These S/N's assume optimal signal processing (via 
matched filters), but it is highly unlikely that optimal processing can be 
achieved. The reasons are two: 

First: Models of galactic nuclei, where massive holes (or other massive 
central bodies) reside, suggest that,iuspiraling stars and small holes typically 



will be in rather eccentric orbits P^""" " I This is because they get injected into 



such orbits via gravitational deflections off other stars, and by the time grav- 
itational radiation reaction becomes the dominant orbital driving force, there 
is not enough inspiral left to strongly circularize their orbits. Such orbital 
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eccentricity will complicate the waveforms and complicate the extraction of 
information from them. In fact, the waveforms are likely to be so complex that 
the number of discrete templates needed in a in one-year coherent-integration 
search for these waves may be enormously larger than modern computers can 
handle. As a result, the signal processing may have to be done coherently over 
short stretches of time with incoherent techniques used to combine successive 
stretches, leading to a significant loss of S/N. Efforts to generate waveforms 
from generic orbits and then scope out this problem are barely getting started, 
as of winter 2001-02. (The waveformii_are fairly well understood for generic 
orbits in the SMBH equatorial plane 1123 and for circular orbits inclined to the 



equator, liHil but inclined generic orbits have not yet been studied, and the 



counting of templates even for equatorial orbits and inclined circular orbits 



has barely begun cl!3. Although the influenc£_ipf radiation reaction on the 



orbital evolution is not yet fully understood e£j, it is probably known well 
enough to scope out the data analysis problem.) 

Second: There may be a significant loss of S/N due to "source confusion" 
— i.e. due to the problem of separating out individual inspiral waveforms from 
the background due to other, weaker inspiral sources. At high frequencies 
this confusion noise could perhaps dominate all other noise sources, just as 
WD/ WD confusion noise dominates at lO""* - 10~^ Hz (Fig. p. Efforts to 
scope this out are also barely getting started. j—. 

The event rates for inspirals into SMBH's are not well known. PhinneyE3 
(relying pia. part on simulations of star-cluster evolution around a SMBH by 
Freitag t2£l) has shown that the detection rate is likely to be dominated by ~ 
lOM© BH's spiraling into ~ lO^M© SMBH's. (Phinney shows that ~ lO^M© 
SMBH's are likely to have larger relaxed stellar cusps around them than either 
smaller or bigger BH's; coincidentally IO^Mq is the optimum SMBH mass for 
LISA. And while inspirals of ~ IOMq BH's probably have a smaller space 
density than inspirals of ^ IAMq NS's, because of their larger S/N, they can 
be detected throughout a volume that is (-\/l0/1.4)'^ ^ 20 times larger.) 

The jestimated inspiral rate for ~ IOMq BH's is ^ 1 — 10/yr out to 



1 Gpc E3'E2il, where the S/N^s would be '--^ 60 to ~ 200 if optimal signal 
processing were possible. If the loss of S/N due to signal-processing problems 
is much larger than an order of magnitude (which it might be), then there 
could be difficulty detecting these inspirals. This motivates exploring a change 
of LISA'S design so its noise curve is lowered by a factor 4. Since such a change 
might significantly increase EISA's cost, it is greatly hoped that a scoping of 
the data analysis problem and a firming of the event rate will reveal that no 
design change is needed! 
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3.5 Collapse of Supermassive Stars 



New and Shapiro E£3 have investigated the possibiUty of substantial GW emis- 
sion from the collapse of supermassive stars (SMS's) to form SMBH's in galac- 
tic nuclei. The formation mechanism for SMBH's is still highly uncertain; one 
path could lead through rotating SMS's. While non-rotating SMS's suffer 
from a well-known radial collapse instability (arising from GR effects), even 



a small amount of rotation can stabilize the star E£3. As the SMS cools and 
contracts, the ratio oi (3 = T/\W\ (ratio of kinetic to gravitational binding en- 
ergy) grows, and (for sufficiently large initial /3) the star eventually becomes 
unstable either to mass-shedding at the equator or to a dynamical, m = 2 
(bar-mode) instability. New and Shapiro show that if viscosity and magnetic 
fields are insufficient to enforce uniform rotation (so that angular momentum 
is conserved on cylinders during the contraction) , then the bar- mode instabil- 
ity is likely to occur first. 

If the bar-mode perturbation grows to order unity and is long-lived, a 
substantial fraction of the star's rest mass-energy could be radiated in GW's. 
These waves end up in the LISA band for a range of plausible initial param- 
eters. For example, for SMS mass M ~ 10^ Mq, initiai4pre-collapse) radius 



R ^ lO^'^cm and initial /3 ~ 10"^, New and Shapiro EH estimate the SMS 
could radiate away ~ 10'^ M© in GW's at / 3 x in a "burst" lasting 

f ~ lyr, which is n = /f ~ 10^ cycles. // optimal signal processing could be 
used to search for these waves, the signal to noise ratio would be S/N ~ hc/hn 
where he = h^/n ^ 10~^^ at IGpc distance; here h is the wave amplitude as 
estimated by New and Shapiro (corrected by them from a different value given 
in their paper), and he is called the waves' characteristic amplitude^^^. Com- 
paring with LISA'S noise curve (Fig. H) we see that S/N ~ 40 for optimal 
signal processing. With sufficiently accurate templates from numerical simu- 
lations of the collapse, it might be possible to achieve a search within a factor 
4 of optimal, so S/N ^ 10, which would be adequate for confident detec- 
tion; but as yet it is far from clear one can do so well — especially as the 
search must be performed in the thick of the WD- WD background. To get 
a clearer handle on the prospects will require considerable analysis — both 
source simulations and data-analysis explorations. 

3.6 LISA and the GW Stochastic Background 

While ground-based detectors can search for a stochastic background by cross- 
correlating their outputs, this method is not available to LISA, since only one 
space-based detector is currently planned. 

However, it is possible to combine the signals from the 3 spacecraft to 
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construct a "symmetrized Sagnac observable" , which is practically insensitive 
to both laser phase noise and gravitational waves Ea. With the aid of this ob- 
servable, one can directly measure (symmetrized averages of) the proof-mass 
noise and optical-path noise, and thereby determine the instrumental noise 
power in the (sum of the) 3 Michelson channels used to detect GW's. Any ex- 
cess "noise" above this power is due to GW's. This method is most effective 
in frequency bands where either the proof-mass noise or optical-path noise 
greatly dominates over the other. With this approach, in a one-year integra- 

ilr" in either the IQ-'^-^ < / < IQ-^-^Hz 



tion LISA could detect ^loo^ 



lO^^ ** < / < lO^i '^Hz bands ti3. Unfortunately, LISA's sensitivity to a 



GW background of early-universe origin is about an order-of-magnitude worse 
than this, since any primordial source with fl{f) <^ 10^^*^ will be "covered up" 
by (the isotropic portion of) the galactic and extragalactic WD/ WD binary 
backgrounds. LISA's low- frequency (LF) sensitivity of ~ 10~^° is about 
two orders of magnitude better than the adjacent-band sensitivities of ad- 
vanced LIGO IFOs (~ 5 X 10"^ in the HF band) and current pulsar timing 



(~ 10"** in the VLF band); Sec. |J 



In Sec. 2.£ we discussed a number of speculative early-universe sources 
that could produce GWs in earth-based detectors' HF band. Variants of 
all these sources could also radiate in LISA's LF band: (i) Universal infla- 
tion (with the standard model predicting waves much too weak for LISA to 
see, n, <^ 10~^^-10~^^, while some less plausible models such as the pre-big- 
bang suggest stronger waves); (ii) Phase transitions of quantum fields (LISA's 
frequency band is ideal for probing the electroweak phase transition, which 
should have occurred at universe age ^ 10^^^ sec and if strongly first order 
could have produced waves strong enough for LISA to see lliil); (iii) Gold- 
stone modes of scalar fields that arise in supersymmetric and string theories; 
(iv) excitations of our universe as a brane in higher dimensions (with LISA 
probing excitation lengthscales ^ 1-10~^ mm if there are one or two relevant 
higher dimensions); and (v) Cosmic strings (with LISA's searches for bursts 
from string cusps and kinks able to constrain string tensions to /i < 10^^'^). 
For further discussion and references, see Sec. 

There is some hope that a follow-on mission to LISA might reach the 
sensitivity required for probing standard inflation, fl{f) ~ 10^^^ — 10^^^. 
To do so will require that the mission's detectors operate outside the regime 
dominated by close binaries. A "short-LISA" with ~ 100 times shorter arms 
would have good sensitivity at / ~ IHz, where the WD binaries disappear 
and the binary background is dominated by NS binaries; a constellation of 3 
short-LISA's spaced evenly around the Sun (along Earth's orbit) could have 
angular resolution A9 ^ 10^^ rad, allowing (probably) the NS binaries to be 
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detected individually and subtracted out. At present, such a constellation of 
9 — 12 satellites seems the best possibility for probing the inflationary GW 
background from space. 



4 Conclusions 

It is now about 40 years since Joseph Weber initiated his pioneering devel- 
opmenL_pf gravitational- wave_detectors IH3, and 30 years since Robert For- 
ward tl3 and Rainer Weiss tl3 initiated work on interferometric detectors. 
Since then, hundreds of talented experimental physicists have struggled to 
improve the sensitivities of gravitational- wave detectors, and hundreds of the- 
orists have explored general relativity's predictions. 

These two parallel efforts are now intimately intertwined and are push- 
ing toward an era in the not distant future, when measured gravitational 
waveforms will be compared with theoretical predictions to learn how many 
and what kinds of relativistic objects really populate our Universe, and how 
these relativistic objects really are structured and really behave when quies- 
cent, when vibrating, and when colliding. By about the time of GR-18, if not 
sooner, we should have some answers. 
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